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ABSTRACT 

The recent discovery of a planetary transit in the star HD 209458, and the 
subsequent highly precise observation of the transit lightcurve with HST, is en- 
couraging to search for any phenomena that might induce small changes in the 
lightcurve. Here we consider the effect of the quadrupole moment of the parent 
star and of a possible second planet perturbing the orbit of the transiting planet. 
Both of these cause a precession of the orbital plane and of the periastron of the 
planet, which result in a long term variation of the duration and the period of 
the transits. For a transiting planet at 0.05 AU, either a quadrupole moment 
similar to that of the Sun or the gravitational tug from an Earth-like planet on 
an orbit of semimajor axis ~ 0.2 AU and a relative inclination near the optimal 
45° would cause a transit duration time derivative of ~ 1 second per year. 

Subject headings: planetary systems - stars: rotation 

1. Introduction 

The first transiting extrasolar planet has been identified around the GO star HD209458 
(Charbonneau et al. 2000; Henry et al. 2000b). The planet has a mass Mp = 0.62Mjup and 
orbits the star with a period of 3.52 days with semimajor axis a = 0.047 AU. Subsequent 
observations of the transit with the STIS camera on board HST have provided highly precise 
photometry of the lightcurve during the transit (Brown et al. 2001), with about 300 flux 
measurements of an accuracy close to 10~^ magnitudes. The shape of the lightcurve has 
been used to determine the radius of the planet, Rp = 1.35Rjup- 



- 2 - 



As discussed by Brown et al. (2001), the determination of transit lightcurves of such 
high precision opens up many opportunities for detecting other phenomena that could induce 
small changes in the lightcurve, such as the presence of satellites or rings around the planet. 
As more transits by other planets are discovered and observed with high precision, the 
opportunities to observe these effects will rapidly increase. This paper considers the effect of 
perturbations of the planet's orbit due to the quadrupole moment of the star or to a second 
planet. These perturbations should cause a precession of the periastron of an eccentric orbit, 
and a precession of the line of nodes if the orbital plane is not coincident with the stellar 
equator, or the plane of the second planet. The precession of the periastron and the line of 
nodes should result in long-term variations of the period and the duration of the transits. 
We will compute these two effects and discuss their detectability with future space missions. 



2. Precession of the orbital plane 

We consider first the effect of the precession of the orbital plane in the simple case 
where the orbit is circular. The orbital angular momentum of the planet is Lp — Mpua?, 
where n = {GMs/a^Y^'^ is the mean motion (or orbital angular frequency), and Mg is the 
mass of the star. We consider two possible causes for the precession of the orbital plane: 
the quadrupole moment of the star, which for solid body rotation has angular momentum 
Lg ~ k'^MgR'^uJs (where Rg is the stellar radius, uj^ the angular frequency of rotation, and 
A;^ ~ 0.1 for main-sequence stars; e.g.. Ford, Rasio, & Sills 1999), or a second planet, with 
angular momentum L2 = M2n2a|, where we will assume 02 > a. The precession occurs 
relative to the plane perpendicular to the total angular momentum vector = Lp -|- L(s 2), 
which will be referred to as the mean plane hereafter. The subindex (s, 2) means that we 
are considering either the stellar rotation or a second planet. 

Let i be the inclination of the orbital plane relative to either the stellar equator or the 
plane of the second planet, and ip and i(^s,2) the inclinations of these two planes relative to 
the mean plane. Obviously, i — ip + i{s,2), and LpSinip — 2) sini(s 2)- We define the 
X-axis as the intersection of the mean plane and the plane perpendicular to the line of sight, 
and as the angle between the mean plane and the line of sight. The line of nodes is the 
intersection of the orbital plane of the planet and the mean plane, and forms an angle Q 
relative to the x-axis. The orbital precession consists of the rotation of the angle Q, with a 
precession angular frequency fl. 
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2.1. Precession induced by the stellar quadrupole moment 

The potential of a star expanded up to its quadrupole moment is given by 



2 



= V— 5 — (Ssm^e*-!), (1) 



r 2 r 

where Rg is the equatorial radius of the star, r is the distance to the center, 6 is the an- 
gle relative to the equatorial plane, and J2 is the quadrupole moment. Estimates of the 
quadrupole moment of the Sun arc quite uncertain and over the range 10^^ to 10^^ (Godier 
& Rozelot 1999; Rozclot, Godier, & Lcfcbvrc 2001). For self-similar main-sequence stars, 
J2 oc UgRl/Mg. In the case of the star HD209458, the spcctroscopically measured rotational 
velocity indicates a rotation slightly faster than the Sun (Queloz et al. 2000). We will use a 
fiducial value J2 = 10~^ in this paper. 

To evaluate the time-averaged torque acting on the planet's orbit as a secular pertur- 
bation, we calculate the potential of a ring of mass Mp in the quadrupole potential of the 
star, 

V = -(3sm ^sm ^ - I) ^ ^-sm z - 1 j . (2) 

The torque is 

3GM,MpRlJ2 ... ... 

T = —dV/di — ^ g — sm2i . (3) 

Since the component of the planet's angular momentum that is changing is LpSmip, the 
precession frequency is Q = t / {LpSmip). We will see below that the quantity determining 
the change in the transit width is Qsinip, which is 

<\ ■ ■ -Rs 3 J2 ... /,s 

\ I smij, — — — n ^r—— sm zi . (4) 
^ Lp A ^ ^ 



2.2. Precession induced by a second planet 

The presence of a second planet will in general cause perturbations on all the orbital 
elements. When averaged over long timcscalcs, the perturbations are classified as secular 
or resonant, and can be calculated with Lagrange's planetary equations (e.g., Murray & 
Dermott 1999, §6). Here we will consider only the secular precession of the fine of nodes, in 
the most simple case when both the perturbing and perturbed planets are on circular orbits, 
which is easily obtained by replacing the planets by uniform rings of mass. Moreover, we 
will use the approximation 02 3> a. The case 02 < a is less interesting, because the most 
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accurate measurements of transits will be for planets very close to their stars; however, the 
precession effects caused by a planet interior to the transiting one can be calculated similarly 
to the exterior case. 

The potential of a ring of radius 02 and mass M2, as a function of the cylindrical radius 
r on the plane of the orbit and the height z above the plane, is, up to second order in r/a2 
and z/a2, 

^ -GM2 / r2 \ , ^ 

a2 V 2a|/ ^ ^ 

The potential energy of interaction of this ring with another ring due to the inner planet, of 
mass Mp and radius a, and an inclination angle i between the two orbital planes, is 

and the torque is 

r = -dV di = sin 2i . (7) 

Just as before, the precession frequency is 

6 • • ^ ^2 a' 3 

sZsmz„ = — — n 5--sm2z. (8) 



2.3. Effect on the transit duration 

The properties of the transit lightcurve are determined by the angle a between the 
orbital plane of the planet and the line of sight. Simple spherical trigonometry shows that 

sin a — sin ip cos (3 cos — cos ip sin (3 , (9) 

where P is the angle between the mean plane and the hne of sight. The duration of the 
transit is 

2{Rs + Rp) , , 

td = — cos 7 , (10 

na 

where the angle 7 is related to the impact parameter b of the transit by {Rs + Rp) sin 7 = 
h — a sin a. The relevant question to determine the observability of the precession of the 
orbital plane is if the rate of change of ta can be measured over a reasonable observing time 
basehne Aio6s ~ 10 years. We will see that the typical value of the precession period is much 
longer than 10 years, so only the time derivative dtd/dt matters: 

dtd a sin 7 da a sin7 • . . ^ . /-.^x 

— td :r = 'td n — smzp cosp smiz . (11) 

dt Rs + Rp COS'' ^ dt Rs + Rp cos"' 7 
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Using ntd — 2{Rs + Rp) 0037/0, we can reexpress this equation as 

dtd Vtsinip . 

— — = ^ 2tan7 cosp smiZ . (12) 

dt n ' ^ ^ 

We can now estimate the value of Cl sin ip and see if the variation of the transit duration 
would be observable. A quadrupole moment of the star J2 ~ 10~^, with a ~ 10ii!s, and 
sin2i ~ 0.1, yields a precession frequency Q sin ip ~ 10~^n, from equation (4). In the case of 
the perturbation by a planet, if its mass is M2 ~ 10~^Ms (i.e., a gas giant), the effect of the 
planet would be comparable to that of the quadrupole for an orbital size 02 — 30a. For a 
51-Peg type perturbed planet with a ~ 0.05 AU, the second Jupiter-mass planet would be at 
a2 ^ 1.5 AU. Such a massive planet should also be detected from the Doppler measurements 
in any case. A more interesting case is if there is a lower mass planet (not detectable by 
Doppler measurements) on a smaller orbit. For example, for an Earth-like planet, with 
M2 ~ 3 X 10~^Ms, at 02 ~ 5a ~ 0.25 AU, and sin2i ~ 0.1, the precession frequency is also 
Qsinip ~ 10~^n. Therefore, from equation (12) (assuming that the trigonometric factors 
are of order unity) , in these cases we should expect a time variation of the transit duration 
dtd/dt ~ 10"^ 

What is the highest accuracy to which we can measure dtd/dt ? In the STIS observations 
of Brown et al. (2001) of HD209458, the flux varies by ~ 1% during a time t^ ~ 0.01 
days at the beginning and the end of the transit. Each one of their flux measurements is 
obtained from a 60 s integration and has a relative accuracy of ~ 10^^. Therefore, each 
data point taken during the falling or rising part of the lightcurvc gives us the starting or 
ending time of the transit to an accuracy ~ (10^^/0. 01)tc, or ~ 10 seconds. Brown et al. 
have measured nearly 100 data points in the rapidly varying part of the lightcurvc, so their 
present data should allow to determine the transit duration to an accuracy approaching 1 
second. Repeating the measurement 3 years later would yield dtd/dt to an accuracy of 10~^. 

The star HD209458 has magnitude V = 7.64, implying that the number of photons 
received by HST over a 60 s integration time is ~ 10^. Therefore, the photometric accuracy 
achieved by Brown et al. could be improved by only a factor ~ 3 by reaching the photon shot 
noise limit; alternatively, a telescope aperture of 0.7 m could reach the same photometric 
accuracy with improvements in the detection efficiency. A dedicated mission observing all 
the ~ 1000 transits over a period of 10 years of stars similar to HD209458 with the same 
accuracy as Brown et al. would yield dtd/dt to an accuracy ~ 10"^'^. 

We can easily change the parameters of the examples discussed previously to consider 
cases where the orbital perturbations would be detectable with the more easily achievable 
accuracy of dtd/dt ~ 10~^. For example, a perturbing planet with Mp = 3 x 10~^Ms, 



-6- 



a2/a — 4:, and sin2i = 0.3, would be detectable at the 5-a level. 



3. Precession of the periastron 



The orbits of the 51-Peg type planets tend to be circular, as expected owing to tidal 
dissipation in the planet, which circularizes the orbit on a timescale (Goldreich & Soter 1966) 



where the factor Q is inversely proportional to the dissipation rate. For the planet Jupiter, 
Q ~ 10^ has been found observationally from the tidal effect on Jupiter's sateUites (loannou 
& Lindzen 1993). Assuming the same Q for HD209458, for which case Mp/M, ~ 5 x 10"^ 
and a/Rp ~ 75, we find Tg ~ 10^ years. How the parameter Q may vary among different 
planets is highly uncertain, because the origin of the dissipation is not well understood (e.g., 
Murray & Dermott 1999, §4.13). So, while internal dissipation in the planets is a plausible 
explanation for the orbital circularity of the closest planets, it is possible that some of the 
closest planets may have eccentric orbits. In addition, the extreme sensitivity of the time 
Te to a/Rp implies that planets that arc only slighly further out from their star can have 
significant eccentricities even if Q is constant. Orbital eccentricities can also be excited 
by planetary companions (e.g., Rivera & Lissauer 2000). In fact, among known extrasolar 
planets, high eccentricities are common when a > 0.1 AU (Butler et al. 2000), and the planet 
around HD217107 has a = 0.072 AU and e = 0.14 ±0.05 (Fischer et al. 1999). We can expect 
that transits of planets in eccentric orbits will be discovered in the near future. 

We note here that tidal dissipation in the star can also circularize the orbit, and make 
the orbit coplanar with the stellar equator. However, the timescale for dissipation in the 
star is much longer than the age of the system, and is in any case similar to the timescale 
for orbital decay (Rasio et al. 1996; Zahn 1977). 



The periastron of an eccentric orbit should precess due to three effects: the relativistic 
precession, and the same two effects considered previously, the stellar quadrupole moment 
and the perturbations from other planets. 

The relativistic precession rate is given by (e.g.. Landau & Lifshitz 1951) 




(13) 



3.1. 



Precession rates 




(14) 



- 7- 



where e is the eccentricity, and we denote the precession of the periastron as uj. For an orbital 
period of ten days, a ~ 0.1 AU, and small eccentricity, this gives ~ 4 x 10~^n ~ 10~'^yr~^. 

For the stellar quadrupole moment, we assume here that the orbital plane coincides 
with the stellar equator. The potential of the star is given by equation (1) with ^ = 0. The 
precession rate for small eccentricities is easily obtained using the epicycle approximation. 
The orbital angular frequency is — {l/r)d(f)/ dr, and the epicycle frequency is (e.g., Binney 
& Tremaine 1987) 

2 dn^ 2 2 SGMsRp2 .... 
K — r— \-4n —n r— ^ — . (15) 



The precession frequency is 



w = n- n , (16) 

where we have substituted a = r and have used zu <^ n. Typically, J2 — 10~^ and Rs ^ 0.1a, 
so UJ < lO^^n. Except in rapidly rotating stars, the effect of the stellar quadrupole is always 
much less than the relativistic precession. 

For the effect from a second planet, repetition of the same method with the potential 
in equation (5) at z—0 yields a precession frequency 

For an Earth-like planet with M^jM^ ~ 3 x 10~^, and 02 — 2a, the precession rate is 
ri7 ~ 3 X lO^^n, comparable to that from the relativistic precession. All the three effects we 
have discussed cause an advance of the periastron, and therefore the total precession is the 
sum of the three effects. 



3.2. Detectability of the periastron precession 

The periastron precession causes both the period and the duration of the transits to 
change. We discuss first the variation of the period. 

In the epicycle approximation (valid for small eccentricities), the eccentric orbit of the 
planet can be described by the motion along an epicycle relative to a circular orbit, with 
coordinates x(t) along the inward radial direction and y{t) in the backward tangential di- 
rection along the orbit. For an orbit in a quasi-Keplerian potential (where the precession of 
the periastron is small), we have (Binney & Tremaine 1987): 



x{f) — aecos{Kt + ipo) 



(18) 
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y{t) = -2aesm{Kt + tpo) ■ (19) 

The central time of the transits is determined by y{t) at the time of the transit. Between two 
successive transits, the epicycle phase — Kt + ipo changes by A'^ = 2Tr{w/n), and therefore 
y{t) changes by Ay — 2ae2Tr{w/n) cos(«;i + ■0o)- Thus, the observed transit period, Pj, is 
(to first order in e) 

A first possible method to detect the periastron precession is to compare the transit 
period Pt with the true orbital period P = 27r/n determined from the Doppler measurements. 
The accuracy of the Doppler measurements are at present eo ~ 0.1 times the velocity 
variation amplitude, and they are unlikely to improve very much due to limitations set 
by photospheric turbulence. Hence, the orbital period can be determined to an accuracy 
AP/P ~ eD-P/(27rA^^/^io), where N is the number of Doppler measurements and to is the 
total time of observation. For = 0.1, N = 1000 and P/to = 10~^, the accuracy achieved 
is AP/P ~ 10~^, allowing detection only for a planet with M2/M5 > 10~^ if e ~ 0.2 and 
02/0 ~ 2. 

A second possibihty is to detect the change of Pt with time: 



dPt . ( VD 



2 



If the time of each transit can be measured to an accuracy of 1 second, the accuracy of the 
period after observing N — 1000 transits would be Is/N^/'^ ~ 10~^-^ s, and for an observing 
time of 10 years the accuracy of dPt/dt would reach ~ 10~^^. For e ~ 0.1, this allows for a 
10-(7 detection of a planet with M2/M5 = 10"^ and 02 = 2a. 

Next we evaluate the change in the duration of the transit, which is caused by the 
variation of the planet's velocity during transit as the periastron precesses. From equation 
(19), the velocity of the planet at transit is v — na{l — 2ecos'0) (notice that for the purpose 
of computing this velocity we can use the approximation k c:^ n), so the duration is 

2{R, + Rp) cos 7 ^ {Rs + Rp) 

td — ~2^ ^ C0S7 (1 + 2ecos'0) , (22) 

V na 

and its time derivative is 

dtd w{Rs + Rp) . 

— = 4e — cos 7 sm . (23) 

dt n a 

As mentioned previously, the perturbation by an Earth-like planet with 02 = 2a implies a 
periastron precession rate zu/n ~ 3 x 10~^; for Rg/a ~ 0.05, 4e ~ 1, we infer dta/dt ~ 10"^, 
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which is not difficult to detect as discussed in §2.3. Prom equations (8) and (12), we see that 
the transit duration change due to the precession of the orbital plane would be comparable 
if sin 2% ~ 0.1. 

In principle, the change in transit duration due to the precession of the periastron and 
the fines of nodes can be distinguisfied by tfie effect on tfie sfiape of tfie ligfitcurve. Precession 
of tfie periastron cfianges only tfie velocity, so tfie sfiape of tfie ligfitcurve does not change. 
However, precession of tfie line of nodes cfianges tfie impact parameter of tfie transit. Of 
course, separating tfie two effects requires a fiigfier pfiotometric accuracy, and will increase 
tfie measurement errors unless tfie effect of tfie periastron precession is also measured from 
tfie cfiange in tfie transit period. Tfie most favorable case would be wfien tfie planet transits 
over tfie edge of tfie star, wfien tfie figfitcurve sfiape fias tfie fastest cfiange witfi a variation 
of tfie impact parameter. 

It is simple to see wfiicfi one of tfie two measurements, tfie cfiange in tfie period or tfie 
duration of tfie transit, can yield tfie fiigfiest accuracy to measure tfie ratio wjn. Tfie period 
derivative is always measured to an accuracy better tfian tfie transit duration derivative by 
a factor P/^o, wfiere is tfie time of observation. Equating tfiis to tfie ratio of equations 
(21) and (23), we find tfiat tfie transit duration cfiange yields a fiigfier accuracy until tfie 
time of observation is to < •ri7~^2(i?s + i?p) cos7/a, and tfie period cfiange gives tfie most 
accurate determination after tfiat. For a close-in planet witfi 2Rs/a ~ 0.1, and an Eartfi- 
mass perturbing planet witfi 02/0. ~ 2, we find zu^^ ~ 3 x lO'' years. Tfierefore, tfie transit 
duration cfiange sfiould yield tfie periastron precession rate witfi tfie fiigfiest accuracy, except 
for substantially more massive perturbing planets. 



4. Conclusions 

Accurate pfiotometry of transit ligfitcurves of extrasolar planets allows tfie detection 
of slow orbital perturbations wfiicfi affect tfie ligfitcurve. Here, we fiave considered tfie 
perturbation due to a stellar quadrupole moment or to a second planet as possible causes for 
tfie precession of tfie line of nodes and of tfie periastron of tfie orbit. Botfi types of precession 
cause a time variation of tfie duration of tfie transit, and tfie periastron precession also causes 
a variation of tfie period. Tfie rate of tfiis variation is given by equations (12), (23), and 
(21). 

Tfie quadrupole moment of tfie star is important for tfie close-in planets, at a ~ 10i?s. 
Tfiese planets are likely to be on circular orbits and tfierefore tfie precession of tfie line of 
nodes sfiould usually be tfie most important effect. Prom observations of tfie stellar spectrum 



-10- 



one can infer the projected rotational velocity of the star, and monitoring the Call flux can 
also reveal the rotation period (Henry et al. 2000a). In addition, spectroscopic variations 
of the star during the transits can in principle be used to measure the angle Q of the line 
of nodes (Quclloz ct al. 2000). If all these quantities arc measured, it should be possible to 
infer from the precession rate of the orbital plane, the quadrupole moment of the star J2, 
which is sensitive to the state of rotation of the inner parts of the star. 

The precession induced by other planets will be added to that caused by the stellar 
quadrupole, and only the sum of both effects is observable. However, once the stellar rotation 
is measured the effect of the quadrupole is predictable to some extent, and it should be 
negligible compared to the effect of another planet in many cases (for large a or M2). 

The precession induced by a second planet depends only on the quantities {M2/ Ms){a/ ., 
and sin 2? (see eqs. [8] and [17]). In principle, by measuring both the orbital and periastron 
precession from the change in the duration and the shape of the lightcurve, or from measur- 
ing the change in the period as well, both of these quantities could be determined assuming 
that there is only one perturbing planet. Even then, the mass and semimajor axis of the 
perturbing planet cannot be separately inferred, and the detected precession could be the 
result of a low-mass planet close to the transiting planet, or a more massive planet further 
away, or a disk of planetessimals. 

In the examples mentioned previously in this paper, we have considered perturbations by 
planets of low enough mass that they would not be detected via the Doppler measurements 
of the star. However, transits could also allow the detection of perturbations by other known 
planets. Several cases of systems with two or three planets have now been discovered, some 
of which seem to be on resonant orbits (Marcy et al. 2001). Measuring the orbital and 
periastron precessions in one of these systems, if they were to show transits, could allow 
the determination of their relative orbital inclination. For a case with M2/MS ~ 10^^, 
02/0 = 2^/^, the precession periods may be as short as hundreds of years, making it possible 
that a whole series of transits may be observed over a reasonable time of ~ 10 years as the hne 
of nodes precesses, and enabling very accurate measurements of the orbital perturbations. A 
more general analysis than that presented here would be required to compute the precession 
rates when the two planets are not coplanar, are both on eccentric orbits, and are possibly 
resonant, including also the effect of changes of the eccentricity on the transit period and 
duration. 

I am grateful to David Weinberg and the Ohio State astronomy graduate students 
Khairul Alam, Chris Burke, Julio Chanamc, Frank Delahaye, Dale Fields, Susan Kassin, 
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